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This paper studies the regularity properties of the density of the 
exit measure for super-Brownian motion with (1 + /3)-stable branch¬ 
ing mechanism. It establishes the continuity of the density in dimen¬ 
sion d = 2 and the unboundedness of the density in all other dimen¬ 
sions where the density exists. An alternative description of the exit 
measure and its density is also given via a stochastic integral repre¬ 
sentation. Results are applied to the probabilistic representation of 
nonnegative solutions of the partial differential equation Au = 

1. Introduction and statement of results. This paper is devoted to reg¬ 
ularity results for the density of the exit measure of super-Brownian motion 
with (1 + /3)-stable branching mechanism from a smooth domain of Exit 
measures of superprocesses were introduced by Dynkin in connection with 
applications to partial differential equations (see in particular [5] and [6]). 
Here we use a stochastic integral representation of exit measures to get pre¬ 
cise information on their regularity or irregularity. As an application we pro¬ 
vide a probabilistic representation for all nonnegative solutions of Au = 
in a smooth domain, in the so-called subcritical case where d < 1 + 2/13. 

Let D be a bounded domain of class in {d>2 throughout this 
work). 11 X € D, we write p{x) = dist(x,D'^) for the distance of x to the 
complement of D. We denote by Af^ the space of all finite measures on D, 
which is equipped with the weak topology. If /i € Aip, supp(/r) denotes the 
closed support of p, which is a subset of D, and we set 

= {fj,€ Mp : supp(/x) C D}. 
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The integral of a function cj) with respect to a measure ^ will often be written 
as (/r, (/)). 

Let f5 G (0,1] and let X = {Xt,t >0) be a super-Brownian motion in D 
with (1 + /3)-stable branching mechanism. To be specific, X is a superprocess 
with branching mechanism 'ip{u) = whose underlying spatial motion is 

Brownian motion in killed when it exits D. The process X is a strong 
Markov process with values in whose distribution will be characterized 

in Section 2. If // G we write for the probability measure under which 
X starts from //. In the first two theorems below, we will consider the case 
where the initial value // G (see, however. Remark 3.1). 

As a special case of the martingale problem recalled in Section 2.3, we 
know that for every twice continuously differentiable function (j) on D, with 
compact support contained in D, 

{Xt, </-) = (/i, (/>) + f\xs, i</.) ds + 

Jo 

where Mt{4>) is a martingale under P^. It will be convenient to use the 
notation 

= ^ J^^[o,t]is)4>ix)M{ds,dx). 

Standard arguments then show that the “stochastic integral” 



ID 


f{s, x) M{ds, dx) 


can be defined for a wide class of integrands / (see Section 2.3 and the 
beginning of Section 3). 

Let X^ be the exit measure of X from D. Note that the usual definition 
of X^ involves the associated historical process, which contains more infor¬ 
mation than {Xt,t > 0). Alternatively, one can proceed as in [6] or [8] by 
defining the superprocess as the collection of all exit measures from time- 
space open sets (these include the measures Xt as special cases). The mea¬ 
sure X^ is a random finite measure supported on dD. We prove in Section 2 
that X^ can be obtained via the following approximation, which is of inde¬ 
pendent interest. For every e > 0, set 


De = {xeD:p{x)>e}, F^ = D\Ds, 


and 

noo 

X^idy)=£~‘^ lF,iy)Xtidy)dt. 

Jo 

Then X^ converges weakly to X^ as e tends to 0, in P^-probability (see 
Proposition 2.1). This shows in particular that X^ is a measurable function 
of (W,f>0). 
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It was proved in [1] and [21] that is almost surely absolutely continu¬ 
ous with respect to Lebesgue measure on dD if and only if2<d<l-|-2//3. 
In the case /3 = 1 and d = 2, more can be said: has a continuous density 

(see [15]). In this work we consider the “stable branching” case; that is, from 
now on we concentrate on the case 0 < (3 <1, and we address the question of 
regularity of the density of the exit measure in dimensions 2 < d < 1 + 2/p. 

Our first theorem provides a stochastic integral representation for the exit 
measure and for its density when it exists. 

Theorem 1.1. Let {P£){x,y),x G D,y £ dD) denote the Poisson kernel 
of D, and let a denote Lebesgue measure on dD. Let y G 

(i) For every continuous function (p on dD, P^-a.s., 

POO P 

(1.1) {X^,(p) = {y,PD(p)+ / PD(p{x)M{ds,dx), 

Jo Jd 

where 

Pd<P{x)= f PD{x,y)(p{y)a{dy). 

JdD 

(ii) Suppose that d < 1 + 2/p. Then, for every y G dD, we may define 
under 

(1.2) X^{y)= f Pn{x,y)iJ.{dx) + f f PD{x,y)M{ds,dx). 

Jd Jo Jd 

The mapping y ^ X^ (y) is continuous in {P ^), for any p G[\,lp P), and 
we have X (y) > 0, P^-a.s., for every y G dD. Finally, 

X^{dy) = X^{y)a{dy), P^-a.s. 

To be precise, we should say in the last assertion that we consider a 
measurable modification of the process {X {y),y G dD). 

We now come to the main result of the present work, which deals with 
the regularity properties of X . For any measurable function /: dD i—> M, 
let ||/||_B denote the essential supremum (with respect to Lebesgue measure 
on dD) of / on the relative open set B G1 dD. 

Theorem 1.2 (Regularity and irregularity of density). Let y G JPlpc- 

(a) If d = 2, the process {X^{x),xGdD) has a continuous modification 
under P^. 

(b) Suppose that 3 < d < 1 + 2/p. Then 

ll^^(■)llu = oo whenever X^(U) > 0, for any open set U C dD, P^-a.s. 
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Obviously, the second part of the theorem remains valid if we replace 
by any version of the Radon-Nikodym derivative of with respect to a. 
Thus, when 3 < d < 1 + 2/P, there exists no continuous density of the exit 
measure. 

The main motivation for studying exit measures comes from their connec¬ 
tions with partial differential equations. A basic result of Dynkin [5] shows 
that the exit measure yields a probabilistic solution of the nonlinear Dirich- 
let problem associated with ^Att = To be specific, for any nonnegative 
continuous function p on dD, the function 

(1.3) u(x) =-logE5^[e“^^^’‘^^], xeD, 


is the unique nonnegative solution to the following boundary value problem 
in D: 


(1.4) 


\Xv = in D, 

V = 4> on dD. 


A major problem is to extend this probabilistic representation to all non¬ 
negative solutions of |Au = in D, and to see that this representation 
induces a one-to-one correspondence between solutions and their traces on 
the boundary (defined in a proper way). This problem was solved in [15] 
in the particular case /? = 1, d = 2. Later, Marcus and VCTon [17] gener¬ 
alized the results of [15] by showing that in the so-called subcritical case 
d< 1-I-2//3, there is a one-to-one correspondence between nonnegative solu¬ 
tions and admissible traces. The next theorem gives a probabilistic formula 
for this correspondence. In order to be able to use the results of [17], we 
restrict our attention to the case of the unit ball. 

We need one more dehnition. The range 7^ of A is defined as the closure 
of the set 


y supp(Ai). 
t>o 


Theorem 1.3. Suppose that d < 1 + 2/P and that D is the unit ball of 
M'^. Let K be a compact subset of dD, and let u be a Radon measure on 
dD\K. The function 


(1.5) u{x) 


logE^^ 


ll{7enA'=0} 




X G D, 


solves the equation ^Au = in Conversely, if u is any nonnegative 
solution of ^Au = in L), there exists a unique pair (K, u) such that the 
representation formula (1.5) holds. 
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As the proof will show, the pair {K, v) can be interpreted as the trace of 
the solution u (defined analytically in [17]). 

Let us emphasize an important point. To make sense of the probabilistic 
representation stated in Theorem 1.3, it is crucial to have chosen a specified 
version of the Radon-Nikodym density of the exit measure. In dimension 
d = 2, we may of course choose the continuous density (as was done in [15]), 
but Theorem 1.3 shows that in higher dimensions the right choice is to 
consider the process X (y) as defined in Theorem 1.1. 

Remark 1.1. In the present work, we do not discuss the quadratic 
branching case f3 = 1. However, our results also hold in that case. Both The¬ 
orem 1.2(a) and Theorem 1.3 are proved in [15] in the case /? = 1. Further¬ 
more, the reader will easily check that the stochastic integral representation 
of Theorem 1.1 is also valid in that case: M should then be interpreted as 
the usual L^-martingale measure associated with super-Brownian motion. 

As a matter of fact, this stochastic representation can be used to simplify 
the proof of the key technical lemma of [15]. 

Let us record some convenient notation for future use. In general, if F 
is a set of functions, we write F~^ for the set of all nonnegative functions 
in F. We use c or C to denote a positive, finite constant whose value may 
vary from place to place. A notation of the form c(a, b,...) means that 
this constant depends on parameters a,b,.... If E is a metric space, let 
B{E) be the corresponding Borel cr-algebra [B{E) will also serve as the set 
of Borel measurable functions on E], We denote by C{E) the space of all 
continuous functions on E and by Cb(Li) [resp. B\,{E)] the space of bounded 
functions in C{E) [resp. in B{E)]. We also denote by Cq{D) the set of all twice 
continuously differentiable functions on D with compact support contained 
in D. Finally, if x G and r > 0, B{x,r) stands for the open ball of radius 
r centered at x. 

The paper is organized as follows. Section 2 recalls basic facts about 
super-Brownian motion and states some preliminary results. Theorem 1.1 
is proved in Section 3, Theorem 1.2(a) is proved in Section 4, and part (b) 
of Theorem 1.2 is proved in Section 5. Connections with partial differential 
equations are discussed in Section 6. The Appendix gives the proof of a 
technical auxiliary lemma. 

2. Preliminaries. 

2.1. Estimates for the Green function and the Poisson kernel. Let {GD{x,y)',x,y G 
D) be the Green function of D and recall that {Pd{x, z); x £ D,z £ dD) de¬ 
notes its Poisson kernel. The functions Go and Pd are continuous on D x D 
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and D x dD, respectively, and they have the following probabilistic inter¬ 
pretation. Let > 0; lixiX S D) denote Brownian motion killed at its 
first exit time from D, and let C, be the lifetime of this process. Then, for 
any 4> £ and x G D, 






(t>{y)GD{x,y)dy, 


nx[<?!>(^C-)] = / 4>{y)PD{x,y)cr{dy). 

JdD 

We will use the following estimates. For every x,y a D and z G dD, 


(2.1) GD{x,y)<C{D)p{y)\x-y\^ 

(2.2) GD{x,y) < G{D)p{x)p{y)\x - 


and 

(2.3) Pd{x,z) <G[D)p{x)\x — z 


Estimates (2.1) and (2.2) can be found in Theorem 2.3 of [22] in dimension 
d>3. In dimension d = 2, they both follow from the more precise bound in 
Theorem 6.23 of [2]. Finally, (2.3) is a consequence of (2.2) and the inter¬ 
pretation of the Poisson kernel as half the normal derivative of the Green 
function at the boundary (see Proposition 5.13 in [2]). 


2.2. Super-Brownian motion and its exit measure. In this section we re¬ 
call the basic facts about super-Brownian motion that will be used in the 
proofs of our results, and we also discuss properties of the associated martin¬ 
gale measure. Without additional effort, the results of this section are valid 
in a more general setting than in the Introduction, namely, for a branching 
mechanism function ^ of the type 

V’(u) = J n{dr){e~'^^ — 1-\-ur), u>0, 

where n{dr) is a cj-finite measure on (0, oo) such that f (r A r^)n(dr) < oo. 
Note that 'il’(u') > 0 for every u > 0. 

Our super-Brownian motion X with branching mechanism ip is a. time- 
homogeneous Markov process in whose semigroup is characterized as 
follows: For every p, G M^,4> G {D) and t > 0, 

=exp-(/r,ut), 


where the function {ut{x),t >0,x € D) is the unique nonnegative solution 
of the integral equation 


ut{x) + Ha; 



'lp{ut-s{^s)) ds 


na;[</’(Ct)l{f<^}] 
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(see, e.g., [6] or Chapter 2 of [16]). In particular, taking (j) = X> 0, we get an 
expression for the Laplace transform of {Xt,l), from which one easily sees 
that E^[(Xi, 1)] < (//, 1) for every t > 0. 

From the preceding Laplace functional, it is not hard to derive that for 
any 








exp-(/r,u). 


where the function (v(x),x £ D) is nonnegative and solves the integral equa¬ 
tion 


v{x) -hlla; 






In view of approximating the exit measure , we now write the following 
joint Laplace transform. For any g G B^{dD) and (j) G B^{D), 


E,. 


exp 


{Xt,<p)dt-{X^,g) 


= ex.p-{g.,w), 


where the function {w{x),x G D) is nonnegative and solves the integral equa¬ 
tion 


W{X) + Hj; 


ijj{w{Ct))dt 


= n. 


dt + g{^(;-] 


This statement is a special case of Theorem 1.1.8 in [6]. We can now prove 
the approximation of the exit measure stated in the Introduction. 


Proposition 2.1. Let X^ he defined as in Section 1. Then X^ con¬ 
verges weakly to X^ as e | 0, in E^-prohahility. 

Proof. Let ip G C^{D). It is enough to prove that 

{X^,ip)^{X^,ip) 

in P^-probability, as e —> 0. To this end, we need only check that, for every 
A,A'>0, 

E^[exp(-A(Xf - A'(X^,(/p))] ^E^[exp(-(A + A')(X^,V^))] 

as e —> 0. We hx A and X' and establish the preceding limit. 

By our dehnition of XP, and results recalled before the statement of the 
proposition, we have 

E^[exp(-A(Xf ,y?) - X'{X^,ip))]=exp{-{ia,wP), 
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where 

(2.4) 

Similarly, 


w^{x) +n,. 


= n. 


K 
Uo 


dt 


= h^{x). 


where 

(2.5) rt;(x)+n. 


E^[exp(-(A + M){X^,ip))] = exp(-(//,'u;)), 
C 


i'{w{it))dt 


= Ii^[{\ + \')if{i^_)\=h{x). 


By standard arguments (see, e.g., the proof of Theorem 1.1 in [5]), (2.5) is 
equivalent to the boundary value problem 

= 'ii){w) in D, 

w = {\ + X')ip on dD. 


Uniqueness of the nonnegative solution for this boundary value problem is a 
consequence of the maximum principle, and so we see that w is the unique 
nonnegative solution of (2.5). 

An application of the bounds (2.1) and (2.2) shows that there exists a 
constant C{D) such that for every x G D and e G (0,1], 




uo 


'^Fe{Ct)dt 


= / GD{x,y)dy<C{D)e^ 


To get this, first note that by the strong Markov property it is enough to 
consider the case when x G F^, and then use the bound (2.1) when \y — x\ <e 
and the bound (2.2) when \y — x\ > e. The point is to observe that the 
Lebesgue measure of Fg r\B{x,S) is bounded above by C\D)e6'^~^ for every 
6 G [e, oo). 

It follows from the previous bound that the functions , e G (0,1], are 
uniformly bounded over D, and by (2.4) the same holds for the functions 
w'^, s G (0,1]. We have then 




tFACtM^t)dt 



GD{x,y)ip{y)dy. 


Using either of the bounds (2.1) or (2.2), and the fact that Pd{x,z) is half 
the normal derivative of the mapping y ^ GD{x,y) at z [in other words, 
GD{x,y) 2p{y)PFi{x, z) when y tends to z along the normal to dD at z], 
we easily get 


lim ^ 

£^o e2 , 


GD(x,y)(p{y)dy = / PD{x,z)p{z)a{dz) 
JdD 
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It follows that h^{x) —> h{x) as e ^ 0, for every x € D. 

Let if be a compact subset of D, and let eo € (0,1] such that p{x) > Eq 
for every x € iL. Denote by Co the first exit time from From (2.4) and 
the strong Markov property at time Co) we get that for every x £ and 
e € (0,eo]) 


w^{x) + 113 ; 



ipiw^iCs))ds 


hl{x), 


where the functions Lq are harmonic on and nniformly bounded. As pre¬ 
viously, this integral equation implies that solves = 'tp{w^) in 

and since the functions are uniformly bounded on D, standard analytic 
arguments (see, e.g.. Theorem 3.9 in [12]) show that the functions are 
equicontinuous on K. At least along a subsequence, we may therefore assume 
that converges to a limiting function w, uniformly on every compact sub¬ 
set of D. By passing to the limit in (2.4), we see that w solves (2.5) and thus 
w = w. We conclude that converges to w, which completes the proof. □ 


2.3. The associated martingale measure. For the results of this section, 
it is convenient to equip the underlying probability space D with the fil¬ 
tration {iFt) generated by X, which is completed as usual with the class of 
.Foo-measurable sets which are P^-negligible for every /r G Aip. All martin¬ 
gales or local martingales will be relative to the filtration (JFt). We will use 
the standard notation AXg = Xg — Xg- for the jump of X at time s (no 
confusion should arise from the use of A also for the Laplacian). 

We hrst recall from [3], Section 6.1 or [10] that X satisfies the following 
martingale problem. For every ip € Cq(D) and every / G C^(K), 

f((Xt,cp}) - f{{Xo,ip)) - i [' f'{{Xg,ip)){Xg,Ag,)ds 

Jo 

(f / + r‘p{x)) 

J 0 D J ( 0 ,oo) 

- /((Ais, if)) - f\{Xg,ip))rip{x))n{dr)Xg{dx)j ds 

is a local martingale. 

From this martingale problem, one easily infers that the jumps of X must 
be of the following type. If s > 0 is a jump time of A, then AA^ = rd^ for 
some r > 0 and x G M'^. More precisely, if J denotes the set of all jump times 
of A, the compensator of the random measure 

•= X! ^(sAXs) 

S^J 
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is given by the following formula. For any nonnegative predictable function 
F on M+ X nxM^, 


( 2 . 6 ) 




J2F{s,u;,AXs 


= E,, 


F{s,uj, ^)N{ds, dfj.) 


where N is the random measure on M_|_ x Aip dehned by 

J G{s, fi)N{ds,dfj,) = j ds J n{dr) J Xs{dx)G{s,rSx)- 

See Theoreme 7 in [10] or [3], page 111. 

Let F be a measurable function on M_|_ x Xip such that for every t > 0, 

\ V 2 n 


(2.7) 


E„ I ^ F{s,AX,f\ <oo. 

^s£jn[o,t] 

Following [14], Section II. Id, we can then dehne the stochastic integral of F 
with respect to the compensated measure N — N, 

[ F{s,i^){N - N){ds,dn), 

JO 

as the unique purely discontinuous martingale (vanishing at time 0) whose 
jumps are indistinguishable from the process lj(s)F(s, AXg). 

We shall be interested in the special case where F(s,/r) = F^{s,^) = 
/ (f){s, x)^i{dx) for some measurable function (p on M_|_ x D. [Some conven¬ 
tion is needed when / \(p{s, x)\^{dx) = oo, but this will be irrelevant in what 
follows.] If (p is bounded, then it is easy to see that condition (2.7) holds. 
Indeed, we can bound separately 

Em 

_ \ S<t 


^{AXs, 1)^1{(AN,,1)<1} j 

^{AXs, l)^l{(AXs,l)<l} 


<E^ 


-I 1/2 


s<t 


'( 0 , 1 ] 


r^n((ir)Ep 


Uo 


{Xs,l) ds 


1/2 


< oo. 


and, using the simple inequality of < (ai + Un)^ for any 

nonnegative reals oi,..., Un,, 

\ 1 / 2-1 


E„ 




<E^ 




1 » 1 } 


. s<t 


'(l,oo) 


rn((ir)Ep 


Uo 


{Xs,l)ds 


< oo. 
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In both cases, we have used (2.6) and the fact that 1)] < {fi, 1). 

To simplify notation, we write 

Mt{(j))= [ [ (l){s,x)M{ds,dx) = [ F^{s,fi){N - N){ds,dfj.), 

Jo Jd Jo 

whenever (2.7) holds for F = F^. This is consistent with the notation of the 
Introduction. Indeed, if 4>{s,x) = ¥?(x) where ip G Cq(D), then by the very 
definition, Mt{(p) is a purely discontinuous martingale with the same jumps 
as the process {Xt , ip). Since the same holds for the process 

Mtiif) := {Xt,ip) - {Xo, ip) - \ f\xs,Aip) ds 

Jo 

(see Theoreme 7 in [10]), we get that Mt{(f)) = Mt{ip). 


3. The stochastic integral representation. We return to the special case 
where 'ip{u) = and thus 


n{dr) 


m+^) 

r(i-/3) 


r ^ ^ dr 


for some (3 G (0,1). 

In this section and in the next two, we fix the initial measure p of our 
super-Brownian motion, and we assume that p G Xip^.. To simplify notation, 
we write P instead of P^ and E instead of E^. 

We need to introduce some notation. Let {pj^{x,y), t > 0,x,y ^ D} be 
the transition density of Brownian motion killed on its exit from D, and let 
{SP, t > 0} be the corresponding semigroup. For any measure v G Xip set 



Pt{x,y)v{dx), 


y € D, t> 0. 


Recall that E[(W,i;^)] = Id every t > 0 and 4> ^ 

(this first-moment formula is easy from the Laplace functional of (Xt , cj)) 
recalled in Section 2). 

For any p > 1, we define the Banach space 

jp ^ X D, SP fi{x) ds dx) 


of equivalent classes of measurable functions with finite norms 


V — 


\f{s,x)I’SPp{x) dxds ) 


\ i/p 


10 JD 

Note that if / does not depend on the “time” parameter s, then 

\ i/p 

\f{x)\PGDKx)dx] , 


'D 
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where 

Goi'iy)^ GD{x,y)v{dx), yeD,ueMp. 

JD 


Lemma 3.1. Let for some p^{l +(3,2). Then the martingale 

Mt{4>) = [ [ 4>{s,x)M{ds,dx), t>0, 

Jo Jd 

is well defined, and bounded in L'?(P) for every q G (1,1+ /3). More precisely, 
for every g G (1,1 + (3), 


<c{f3,p,q){U\\l+U\\l). 

Moreover, for any sequence of functions {(pmn > 1} such that (/>„ —> (/> in 
as n—> oo, we have 


(3.1) 


E 


sup|Mi((/))|'' 

.t>Q 


(3.2) 


lim E 

n^oo 


suplMficpn) 

. t>o 


= 0 


Vg G (1,1 + (3). 


Proof. To see that the martingale Mt{4>) is well defined, we need to 
verify condition (2.7) with F = F^. We will in fact prove more by checking 
that, for every g G (1,1 + (3), 


(3.3) 


E 


Y,F{s,AXs)^ 


9/2-1 


< OO. 


First note that since p/2 < 1, we have A whenever 

flj > 0 for every i G /. We use this in the second inequality below: 

AsgJ / 


<E 


<E 


= E 


VseJ 

E1{(AW,i><i}A(s,AX,)P 

S^J 


p/2- 


q/p 


g/p 


'( 0 , 1 ] 


ds J n{dr) j Xs{dx)l^j.<i}X^\4’{s:x)\^ 
r^n{dr 


-\g/p 


r°° f n ^ 

I ds / dx SYy{x)\(l){s,x)\^ 

0 Jd 


= C{/3,p,q)\\cl)\\ 


9 

P’ 
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using (2.6) and the fact that /^q rPn{dr) < oo since p> 1 + p. 
Similarly, 


E 


^{{AXs,i)>i}F'{s,AXs 


q/2' 


<E 

= E 


Y.l{iAX.,i)>i}\F{s,AX,)P 

-S^J 


ds / n{dr) / Xs{dx)l^^yi^r'^\(l){s^x)\^ 


r'^n{dr] 


ds I S^p{dx)\cj){s,x)P 


ID 


'(l,oo) 

= C{P,qm\l, 

using (2.6) and the fact that r^n{dr) < oo since q<l + p. 

By combining the last two bounds, we see that (3.3) holds. Furthermore, 
by the Burkholder-Davis-Gundy inequality for purely discontinuous mar¬ 
tingales (see, e.g.. Chapter VII of [4]), 


E 


sup \Mt{p)P 
Lt>0 


< C{q)E 


J2F{s,AXs 


ql2- 


\ 5Gc/ 


and the bound (3.1) follows from the previous inequalities. The last assertion 
is immediate from (3.1), observing that (/>„—>(/> in implies P 

since the measure p^x) dx ds is hnite. □ 

The next lemma is a Fubini-like theorem for our stochastic integrals. 

Lemma 3.2. Let {E,S,h') be a a-finite measure space and let p he a 
measurable function on M_|_ x D x E. Assume that for some p G (1 + P,2), 

r roo r 

/ / / \Pis,x,y)\PS^p{x)dxdsi'{dy) <oo, 

J E J 0 J D 

and for every y € E, 

POO r 

/ / \p{s,x,y)\^S^p,{x) dxds < oo. 

Jo Jd 

For every y € E set py{t, x) = p{t, x, y) and 

Mt{py)=[ [ py{s,x)M{ds,dx). 

Jo Jd 
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Then, for every t G [0,oo], the process {Mt{(py),y G E) has a measurable 
modification, and 

(3.4) J^Mt{4>y)u{dy) = J^(^J^(j){s,x,y)n{dy)^M{ds,dx), P-a.s. 

Proof. We only sketch the arguments. First note that our integrability 
assumptions guarantee that the stochastic integrals Mfifiy) are well dehned 
for every y € E, that the function y ^{s, x, y) is i/-integrable pL{x) dx ds- 
a.e., and that the stochastic integral in the right-hand side of (3.4) is well 
defined, independently of the value we give to f^(j){s,x,y)u{dy) when y—> 
(p{s,x,y) is not z^-integrable. By standard arguments, it suffices to prove 
the lemma when is a finite measure and = 1a is an indicator function 
(note that the integrability assumptions of the lemma are then automatically 
satished). In the particular case where 74 = ^41 x A 2 , with Ai G .B(M+ x D) 
and A 2 € £, the various assertions of the lemma are immediately verified. 
The general case follows from a classical monotone class argument. □ 

Proof of Theorem l.l(i). Let & C{dD). We may extend i;^ to a 
continuous function on D, which we still denote by cj). By standard tech¬ 
niques (see Proposition 2.13 in [11] or Exercise II.5.2 in [19] for the finite 
variance branching case), it is easy to obtain that for every t>0, P-a.s., 

(3.5) = f f S^_sfi{x)M{ds,dx). 

Jo Jd 

We then apply Lemma 3.2 to the (bounded) function (s, x, t) —s- l^s<tjS[fg4>{x), 
noting that /g°° l^g^^jS[fg(p{x) dt = GD<t>{x). It follows that 

POO POO P 

(3.6) / {Xt,(l))dt= {y,GDfi) + / GDfi{x)M{ds,dx). 

Jo Jo Jd 

From the definition of X^, we get for any e > 0, 

POO P 

(3.7) {X^,fi) = {y,GDn+ GDr{x)M{ds,dx), 

Jo Jd 

where /^(x) = e~‘^lp^{x)fi{x). As in the proof of Proposition 2.1, it is easy 
to verify that, for every x € D, 

(3.8) GDf{x)^PDfi{x) 

as e —> 0, and furthermore, the functions Gpf^ are uniformly bounded 
over D. By dominated convergence, we see that Gpf^ converges to Pnfi 
in for every p G (1 -|- /3,2). By passing to the limit e —> 0 (using the last 
assertion of Lemma 3.1), we get the desired result. □ 
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Proof of Theorem l.l(ii). Let p G (1, ^). From the bounds (2.1) and (2.3), 
it is straightforward to verify that, for any compact subset K of D, 

(3.9) sup (f GD{x,y)PD{y,zf dy) <oo. 

xGK,z£dD\^ D / 


We are assuming d < 1 + 2//3, or equivalently 1 + /3 < We can thus 
choose p G (1 + /3, A 2) and the preceding estimate implies that the 
(time-independent) function {s,y) — > PD{y,z) is in for every 2 : G dD. 
In particular, the stochastic integral appearing in the definition of (z) 
is well defined according to Lemma 3.1. Furthermore, using the fact that 
y G we can apply Lemma 3.2 to the function (y, s, z) Poiy, z) and 

the measurable space {E,£,u) = {dD,B{dD),a). It readily follows that the 
process (X^(z),z G dD) has a measurable modification, and that, for any 
cj)eC{dD), P-a.s., 


{X^,(j)) = {lz,PD(j)) + [ [ ( [ PD{x,z)(p{z)cr{dz)]M{ds,dx) 

Jo Jd \JdD J 


= (j){z) Pd{x, z)pL{dx'^ a{dz) 

-|- [ 4>{z)( [ f PD{x,z)M{ds,dx) ]a{dz) 
JdD \Jo Jd 

= / (j){z)X^ {z)a{dz). 

JdD 


This is enough to conclude that X^{dz) = X^ {z)a{dz), P-a.s. 

In particular, we must have X^(z) > 0, a{dz) a.e., P-a.s. From the esti¬ 
mate (3.9) and the last assertion of Lemma 3.1, it is easy to see that the 
mapping z ^ X (z) is continuous in L'?(P), for every q < 1 + P, and it 
follows that X^{z) > 0, P-a.s., for every 2 ; G dD. □ 


Remark 3.1. The proof of part (i) of Theorem 1.1 does not depend 
on the assumption y G and the result is indeed true for an initial 

measure y G Things go differently for part (ii): If ^ G Al^ \ Al^^,, the 
function (s,x) —> Pd{x,z) may no longer be in for any p > 1 + 0, and 
the stochastic integral appearing in (1.2) may not be defined. Still from the 
additivity property of superprocesses, we can recover from the particular 
case y € the fact that the exit measure is absolutely continuous with 

respect to Lebesgue measure on the boundary. 


4. Continuity of the density in two dimensions. In this section, we as¬ 
sume that d = 2 and we prove part (a) of Theorem 1.2. As we want to use the 
Riemann mapping theorem, we will first assume that D is simply connected. 
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The first term in the right-hand side of (1.2) is obviously continuous in 

y. So, to prove the existence of a continuous modification of X {y), it is 
enough to check the existence of a continuous modification of the stochastic 
integral 

poc r 

Z{y)^ / PD{x,y)M{ds,dx). 

Jo Jd 

Before we continue, let us introduce the following notation. Let Dq be the 
unit disc of the plane and denote by (Jo{dy) the Lebesgue measure on the 
unit circle ODq. The Poisson kernel in this case can be computed explicitly: 

1 1 — IxP 

(4.1) Poix,y) = ^ |^_3.|2 ’ xGDQ,y€dDo. 

The next lemma is crucial for estimating the moments of increments of Z{-). 


Lemma 4.1. 


(4.2) 


(a) Set po{x) = dist{x,dDo). Let a>0, p £ {0,2 + a), and 
2 + a-p, 


7 = 


. r 2 + a 

if —-— <p<2 + a, 


2 + a 


-e, 


ifp = 


2 + 0 


P, 

\ 


if 0 <p < 


2 + 0 
2 


where ££ (0,^^) is arbitrary. Then there exists a constant c = c{p,a,e) 
such that 

(4.3) [ po{xT\Po{x,yi)-Po{x,y2)\^dx<c\yi-y2r Vyi,y2e9T>o- 
Jdo 

(b) For any B C Dq such that dist(il, dDo) > 0, there exists c = c{B) such 
that 

(4.4) sup |Po(a;,yi)-.Po(a;, 1 / 2)1 <c|yi-?/ 2 | yyi,y 2 ^dDo. 

x£B 


The proof of Lemma 4.1 appears in the Appendix. 

Since H is a bounded simply connected domain in M^, the Riemann map¬ 
ping theorem allows us to find a conformal mapping ifi from Dq onto D. 
Under our assumption that D is of class C^, if extends to a one-to-one con¬ 
tinuous mapping from Dq onto D. In fact, we can say more. According to 
Chapter 3 of [20], if' also has a continuous extension to Dq and if' does 
not vanish on Dq. In particular, \if'\ is bounded below and above on Dq 
by positive constants. It is also easy to check that for every x,y £ D and 
z£dD, 

(4.5) 

(4.6) 


PD{x,z) = \lf'{lf ^(z))] ^PQ{lf ^{x),lf ^(z)) 
GD{x,y) = Go{if~^{x),if~^{y)). 
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Let 9{x) = ^( 3 ;))! for every x G D. 

Lemma 4.2. Let p G (0,3), e G (0, |) and p G There exists c = 

c{p, e, D, p) such that for every yi,y 2 G dD, 

[ ivi)P d{ x,yi) -9{y2)PD(.x,y2)fGDp(.x)dx 

JD 

(c\yi-y2\^, i/0<p<|, 

< I c|yi ifp = l, 

[c\yi-y2\^~^, i/|<p<3. 

Proof. Let p be the image of p under and set £{p) = supp(/i), 
= dist{£{p),dDQ), and £{p)f‘*P = {x G DQ:dist{x,£{p)) < p^/2}. Then, 
using (4.5) and (4.6), 

/ \0{yi)PD{x,yi) -9{y2)PD(.x,y2)fGDp(.x)dx 
JD 

= f \Po{x','ilJ~^iyi)) - Po{x',ip-^{y 2 ))f 

JDo 

X GQ{w,x')p{dw)^\'ip'{x')\‘^ dx'. 

By an application of the Fubini theorem, 

/ ([ Go{w,x')p{dw))\'4!'{x')\‘^ dx'^cip,^;), 

JDo \JDo J 

and on the other hand, the bounds (2.1) easily imply that for every x' G 

Do\£{~pY*/\ 

/ Go{w,x')p{dw) < c{p,Y)poix')- 

JDo 

It readily follows that 

/ \0{yi)PD{x,yi) - 9{y2)PD{x,y2WGDp{x)dx 

JD 

<c{p,Y)( sup \Po{x','ilj~^{yi)) - Po{y','ip~^iy 2 )W 

^x'££{P)p*/^ 

+ [ ^ \Po{x',Y~^{yi)) - Po{x',Y~^{y 2 Wpo{x') dx'^ 

JDo\£{py*/^ J 

< c{p,Y){\Y~^{yi) - Y~^{y2)f‘ + \Y~^{yi) - Y~^{y2)V)^ 


18 


J.-F. LE GALL AND L. MYTNIK 


where 7 is as in Lemma 4.1 with a = 1, and we have used both assertions 
of this lemma to derive the last inequality. Since ■0“^ is Lipschitz on D, the 
bounds of the lemma follow easily. □ 


Lemma 4.3. Let p € {1 + (5,2) and q G {1,1 + (5). There exists a constant 
c = c{P,p,q, D, p) such that 


n\e{yi)Z{yi)-e{y2)Z{y2)\^]< 


c\yi -1/2I'', 


ifQ<(5<\, 

if\<l5<l. 


Proof. Recall from the proof of Theorem l.l(ii) that the function 
{s,x) PD{x,y) belongs to for any y € dD and r G (1,3). From Lemma 
3.1, it follows that 

no{yi)Z{yi)-e{y2)z{y2)\‘^] 

< c{P,P,q)(^(^J^ \e{yi)PD{x,yi) - e{y2)PD{x,y2)\^GDh{x) dx'^ 

+ \ 0 {yi)PD{x, yi) - 0{y2)PD{x, y2)\'^GDh{x) dy^ . 

In the case 0 < /3 < ^, choose p G (1 + /3, |) and immediately get the desired 
bound from Lemma 4.2. Similarly, in the case ^ < /3 < 1, the desired result 
follows from Lemma 4.2. □ 


Proof of Theorem 1.2(a). We apply the Kolmogorov criterion of 
continuity to get the existence of a continuous modification of the process 
9{y)Z{y) [and hence also of Z{y)]. The needed bounds for moments of in¬ 
crements of 9{y)Z{y) are obtained from the preceding lemma: In the case 
0 < /3 < ^, this is immediate since q> 1, and in the case | < /3 < 1 , we ob¬ 
serve that we can choose p and q sufficiently close to 1 -|- /? to ensure that 
g(3 — p)/p> 1. The existence of a continuous modification of the process 
Z{y), together with the remarks of the beginning of this section, completes 
the proof of part (a) of Theorem 1.2, in the simply connected case. 

The general case when D is not simply connected can be treated via a 
localization procedure analogous to Section 4 of [15]. Instead of the special 
Markov property of the Brownian snake used in [15], one uses the Markov 
property of superprocesses in the form stated in Theorem 1.1.3 of [ 6 ]. Details 
are left to the reader. □ 


5. Irregularity of the density in high dimensions. In this section, 3 <d < 
1 + 2//9. If z G dD and r > 0, we denote by BQ{z,r) the open ball centered 
at z and with radius r in dD: Bg{z, r) = {y G dD :\y — z\ < r}. 


EXIT MEASURE OF SUPER-BROWNIAN MOTION 


19 


In order to prove part (b) of Theorem 1.2, it is enough to verify that the 
property 

(5.1) ||Al^(-)||b = oo, P-a.s. on the event {X^(B) > 0} 

holds whenever B is a fixed boundary ball. 

We thus fix a boundary ball B = Bg{zo,r]Q). For technical reasons, we also 
introduce a smaller closed ball B' = Bg(zo,??o)> with tjq < rjQ. If SB' denotes 
the relative boundary of B', we assume that cr(SB') = 0 (this is certainly 
true for all but countably many values of t/q). We consider a sequence (en) 
of positive numbers decreasing to 0. For definiteness we may take = 2“”". 
Then, for every integer n > 1, we set 

Bn = {x £ D: dist(x, B') < e„}. 

Lemma 5.1. We have 

noo 

/ Xs(B„) ds —> X^(B') asn—^oo, in F-probability. 

Jo 

Proof. From (1.1), we have for every (p £ Bb('9id), 

= f h{dx) [ a{dy)PD{x,y)ip{y). 

J JdD 

Taking ip = IgWi we see that our assumption cj(SB') = 0 implies X^{dB') = 0, 
a.s. The statement of the lemma is then an easy consequence of the weak 
convergence of towards X^ (Proposition 2.1). □ 

WehxaG(2/(/3 + l),2). Let 

Tn = inf{s > 0: AX^(B„) > e"}. 

Lemma 5.2. We have 

(5.2) P(rn = oo I X^(B') > 0) —> 0 asn—>oo, 

and 

(5.3) limsupP(r„ = oo) <P(X^(B') = 0). 

n—^oo 

Proof. Equation (5.3) is an immediate consequence of (5.2). To verify 
(5.2), we will follow the lines of the proof of Lemma 4.1 of [18]. Define 

Z^ = Ni[0,t]x{p£M^:p{Bn)> £“}), 

where N is the point measure of jumps of the process X, which was intro¬ 
duced in Section 2.3. Then 

(5.4) {r„ = oo} = {Z- =0}. 
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Recall (2.6) for the compensator of N. From a classical time change result 
for counting processes (see, e.g., Theorem 10.33 in [13]), we get that for each 
n there exists a standard Poisson process > 0) such that 

Zf = A- (^c(/3)e-“(^+') ’ 

where c(/3) = /3/r(l — /3) > 0. Fix <5 > 0 such that 2 — a{P + 1) + 5 < 0. Then 

p(z^ = o,a:^(bO >0) 

<p(^A-(e-^) = 0, c(/3)e-“(^+')^'"x,(B„)ds>e-'5, X^(B')>o) 

(5.5) + p(^c(/3)e-“(^+^) Xs(Bn) ds < X^(B') > o) 

< P(A"(e-'5) = 0) 

+ p(^c(/3)4-“(^+i)+' 5 (^e-2 Xs{Bn) ds^ < 1, X^(B') > o) . 

The first term on the right-hand side of (5.5) is P(^"(e“‘^) = 0) = exp 
which converges to 0 as n —> oo. Now, by Lemma 5.1, /q°° Xs(B„) ds 

X^{B'), in probability, as n —> oo. Since 2 — a{j3+l) + 5 < 0, we immediately 
get that 

p(c(/3)4-“(^+i)+^(e-2^°°X,(B„)d5) <1, X^(B')>o) ^0, 
as n —> oo. Hence, the result follows from (5.4) and (5.5). □ 


In order to get a lower bound for X^ in terms of X ^, we observe that there 
exists a positive constant Ci = Ci{D) such that cr(Ba(a;,2e„)) < Cie'^^ for 
every n > 1 and x G dD. If n is large enough so that 2e„ < Vo ~ Vo ^ which 
we assume from now on, we have Bg{z, 2e„) C B for every z G B', and so 



))<CieL‘ 



B- 


Thus 


(5.6) Efexpl-CiljX^llB}] <E 


exp 


■sup 4 '^{X^ABg(z, 2 er.)) 

zGB' 


On the event {xn < oo}, denote by Cn,rn the spatial location and the size of 
the jump at time Tn, meaning that XX^^ = From the strong Markov 

property at time r^, together with the additivity property of superpro¬ 
cesses, we know that conditionally on {r^ < oo}, the process {Xr^j^tA > 0) is 
bounded below in distribution by (Xp, t > 0), where X^ is a super-Brownian 
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motion with initial value rn5c_^. From our approximations of the exit mea¬ 
sure, it follows that conditionally on {xn < oo}, is bounded below in 
distribution by the exit measure of X'^ from D. Hence, from (5.6) we 

get 


Efexpl-CillX^llB}] 

<E 


(5.7) 


l{r„<oo} exp<^ - sup {X , lBg{z,2e^)) 
z&B’ 


+ E(Tn = oo) 


< E 




exp<^ - sup (X , 

' z&' 


+ r{Tn = oo). 


Note that, on the event {r^ < oo}, we have and Cn £ B„. We now 

claim that 


(5.8) lim sup Er 5 


exp<{-supe„ {X ^Bgiz,26^)) 

zSB' 


= 0 . 


To verify (5.8), let xq G and r > e“. By the definition of B„, there exists 
yo G B' such that \yo — xo\ < Sn- Then, using the Laplace functional of the 
exit measure as recalled in Section 2.2, 


Er-5j. 


exp<^-supe„ {X ,lBg{z,2en)) 

' zGB' 


(5.9) 


< Er5^ 


XQ 


exp<{ -ei '^(X^,lBa(j/o, 2 e„)) 


= exp(-r'(;” (xo)) 

<exp(-e>" (xo)), 

where the nonnegative function {Vy^{x),x G D) solves the integral equation 

(5.10) Vy^{x)+ [ GD{x,y)Vy^iy)^^^dy = sl~‘^ f PD{x,z)a{dz). 

JD JBg{yo,2£n) 

Lemma 5.3. Under the conditions 2/(/3 -|- 1) < a < 2 and 3 < d < 1 -|- 
21(5, we have 


(5.11) 


lim 


inf 


oo V3;oGB„,j/oGBg|j/o-a;o|<£i 


e>^o(®o)) =+oo. 
^ / 


Let us postpone the proof of Lemma 5.3. Our claim (5.8) readily follows 
from (5.9) and (5.11). By passing to the limit n ^ oo in the right-hand side 
of (5.7), and then using Lemma 5.2, we arrive at 




E[exp{— ||X'^||b}] < limsupP(rn = oo) < E(X‘^(B') = 0). 
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We can now let B' increase to B by varying along a suitable sequence 
increasing to t/q- Since the event {X^(B) = 0} is the decreasing limit of the 
events {X^(B') = 0 } along this sequence, we get 

E[exp{-||X^||B}] < P(X^(B) = 0). 

Since obviously ||Al^||b = 0 on the event {X^(B) = 0}, the desired property 
(5.1) follows from this last bound. This completes the proof of part (b) of 
Theorem 1.2. 


Proof of Lemma 5.3. Let n > 1 and xq € B„, j/o £ such that 
— yo\ < Gn- In what follows we will need to assume that n is sufficiently 
large, but our bounds will then be uniform in xq and i/q. To simplify notation 
we write u"' = Vy^. Note that by (5.10), for every x € D, 

v'^ix) < f Pd{x, z)a{dz). 

JBg{yo,2en) 

Therefore, 

/ GD{xo,y)v^{yf^^ dy 
, , JD 

(5.12) 

r / r \ 

<ga-rf)(i+/3) / GD{xo,y)[ Pn{y,z)a{dz)] dy. 

We first get a lower bound on the right-hand side of (5.10) for x = xq. 
Since D is of class there is a number a > 0 such that, for every 2 G dD, 
there exists an exterior sphere of radius a tangent to dD at z. Suppose that 
n is large enough so that < a, and for z G dD denote by the closed 
ball with radius Snl’l tangent to dD at z and such that B^ PiD = 0. Then, 
if X G H is such that |x — z| <en, the probability that a Brownian motion 
started at x exits the domain H at a point of Bg{z, 2e„) is bounded below by 
the probability that this Brownian motion hits B^ before exiting B{z,2en)- 
Clearly, this probability is bounded below by a constant Co{d) > 0. Hence, 


(5.13) f PDixQ,z)a{dz) = e BQ{yo,2en))>CQ{d). 

JBg{yo,2en) 


We then turn to an upper bound for the integral over D in the right-hand 
side of (5.12). It will be convenient to deal separately with the integrals over 
D n B(yo, £n) P Id B{yQ, respectively, where 0 < 7 < 1 is chosen so 
that 


d< 


1 -1-7 


+ 1. 
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With obvious modifications, we can then follow the calculations of ([1], 
page 81) and, using (2.1) and (2.3) in the first inequality below, we obtain, 
for n large enough, 


Bg{yo,2£n) 


= [ GD{xo,y)( [ PD{y,z)a{dz)] dy 

J Dr\B{yQ,eZY \JBg{yo,2£n) / 

/ /■ \ 1+/3 

<c{D)[ / aidz) 

\Je ' 

Jl 


\xo-y\^-'^piyf+^ sup \y - z\-'^^^+^Uy 

DnB{yo,£ZY zeBg{yo, 2 £n) 


/ p{yf+f^ sup \y - zl-'^^^+I^Uy 

JDnB(yQ,eZY z£Bg{yo,2£n) 

< / (dist(y, Baivo, 2e4))"+^-''('+^) dy 

JDnB{yQ,elY 

/*diam D 

where the last inequality holds because d < 1 + 2/P implies 1 + /3 — d/3> —1. 

Let us turn to the integral over D fi B{yQ,£/), which is denoted by . 
Notice that, for y G L), 

/ PD{y,z)cr{dz) =ny(^f_ € Ba{yo,2en)) < 1 . 

J Bg{yo,2£n) 

Hence, using again (2.1) and (2.3), 


/-< 


lDnB{yo,£Z) 


GDixo,y) 


< c{D) [ 

JL 


Bg(yo,2en) 
1-d 


PD{y+)(y{dz) dy 


DnB{yo,£Z 


\xo-y\ p{y) 


X 1 


4p(y)<4£n} 


+ lMy)>4£„} / p{y)\y-z\ '^cr{dz)]dy 

JBg(yo, 2 en) J 


<c{D) [ ^ \xo-y\^ '^p(y)(i{p(j/)<4e„} + i{p(y)>4£„}p(y)^ "^4 ^)dy 

JDr\B{yo,£l) 

<c{D)en [ \xo-y\^~‘^dy 

J Dr\B(yo,£Z) 

< c{D)£ne/. 
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By combining the preceding bounds, we get 

/ GD{x,y)v^{yf^Pdy < + ^2 ) 

JD 

(5-14) ^ ^ 1 + 7 ^ 

= c(ei^-'^)^ + e2+/3+7-d-<i/3). 

Therefore, by (5.10), (5.12)-(5.14), we have 

(5.15) v^ixo) > Co(d)4“" - c(Il)(eii-‘')'^ +e4^+^-''-"4. 

Hence, 

(5.16) ey^ixo) > £“+'-"(^^0(4 - c(Zl)(4"-')('-") +4+/?+7-<i/3)) 

for n large enough. Since d < + 1 and 7 < 1, the expression in brackets 

converges to Co((i) > 0 as n —> 00 . Moreover, since d > 3 and a < 2, we have 
—> +00 as n —> 00 , and the desired result follows. □ 

6. The probabilistic representation of solutions of In this 

section, we concentrate on the case when D is the unit ball of M'^, and 
we prove Theorem 1.3. Before starting the proof, let us observe that our 
definition of the range (which agrees with [ 6 ]) is slightly different from the 
one in [7] or [ 8 ]. The reason is that a superprocess is defined in [7] or [ 8 ] as 
the collection of its exit measures from space-time open sets. It is, however, 
not hard to see that both definitions give rise to the same random closed 
set, P^-a.s. for any y G 

We first recall the definition of the trace of a solution following [17]. Let 
u be a nonnegative solution of the partial differential equation 

( 6 . 1 ) = in L). 

We dehne the trace tr(u) of u on the boundary as the pair (iL, v), where K 
is a compact subset of dD and is a Radon measure on dD \ K, which is 
determined as follows: 

(i) A point y G dD belongs to K if and only if, for every relative neigh¬ 
borhood 17 of y in dD, 

lim / u{rz)a{dz) = 00 . 

rnJu 

(ii) For every continuous function y? on dD, with compact support con¬ 
tained in dD \ K, 

lim / u{rz)(p{z)a{dz) = / (p{z)i'{dz). 

rn JdD JdD\K 
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Under the condition d< 1 + 2/(3, Marcus and Veron [17] proved that the 
mapping u —> tr(u) induces a one-to-one correspondence between the set of 
all nonnegative solutions of in D and the set of all pairs {K, u), 

where K is a compact subset of dD and is a Radon measure on dD \ K. 
(In the special case (3 = 1, this result was obtained earlier in [15].) 

Let us prove the first assertion of Theorem 1.3. If u is given by (1.5), we 
aim at proving that u solves (6.1). This is basically a consequence of the 
known connections between superprocesses and partial differential equa¬ 
tions. Consider first the case when i'(dy) = g{y)a{dy), where s' is a non¬ 
negative continuous function on dD, with support contained in dD \ K. 
The random variable Y such that Y = -|-oo on the event {TZ CiK ^ 0} and 
Y = {X^ ,g) on {TZ fl RT = 0} is a stochastic boundary value in the sense of 
[7] (see, in particular. Theorem 6.1 in [7]). Therefore the function 

x=>- loglE^ [exp -Y] = - logE^^ 

solves = in D. 

Coming back to the case of a general Radon measure v on dD \K, we 
may find a sequence of nonnegative continuous functions gn, with support 
contained in dD \ K, such that 

lim / ip{y)gn{y)(j{dy) = / ^p{y)i^{dy) 
n^ooj J 

for every ip € C{D) with compact support contained in dD \ K. 

Lemma 6.1. On the event {TZ r]K = 0}, we have 

{X^,gn)^ J X^{y)v{dy) 

as n=>- oo, in -probability for every xq G D. 

Proof. Let e > 0 and = {y G dD: dist(y, K) < s}. Since 7?. is a closed 
set, the event {TZCiK = 0} is the union of the events {TZCiK^ = 0} over all 
e > 0. Also, on the event {TZ H = 0}, it is easy to see that X^ puts no 
mass on (use Proposition 2.1) and that X^(y) = 0 a.s., for every y G K^. 

Fix e > 0 and let h^:dD [0,1] be a continuous function such that 
hsiy) = 0 il y G K^j 2 and hs{y) = 1 if y ^ iLg. In view of the preceding 
remarks, the proof of the lemma reduces to checking that 

lim {X^,hegn)= [X^{y)he{y)v{dy), 

n^oo J 

in P^^^-probability, for any xq G D,£ > 0. 
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As a special case of ( 1 . 1 ), we have P^^^-a.s., 

PCX) P 

{X^,hegn) = PD{he9n)ixo)+ / PD{hegn){x)M{ds,dx). 

Jo Jd 

Now, for every x € D, we have 

PD{x,y)hs{y)i'{dy), 


PD{hegn){x)= / PD{x,y)he{y)gn{y)(7{dy) 
JdD 


IdD 


as n —> oo. Recall from the proof of Theorem 1.1 (ii) that the (time-independent) 
functions {s,x) Pd{x,z) are bounded in when z varies in dD for any 
p G {1 + (5,{d + l)/{d — 1)). It follows that the previous convergence holds 
in lA for any p G (1 + /3, (d + l)/(d — 1)). By Lemma 3.1, we conclude that 
(A^, /lePri) converges in ), for every <7 G ( 1,1 + /3), toward 


IdD 


PD{xo,y)he{y)v{dy) + 


D \JdD 


PDix,y)hs{y)iz{dy) ]M{dsdx) 


= [ {y)he{y)y{dy), 

JdD 

thanks to (1.2) and the “Fubini theorem” Lemma 3.2. □ 


We come back to the proof of Theorem 1.3. For every n > 1 let 


Un{x) 




ll{7^nA=0} exp 


gn{y)X^{dy)^ , 


X G D. 


We already saw that Un solves (6.1), and by the lemma, Un{x) converges to 
u{x) as n —> 00 , for every x G D. Since the set of nonnegative solutions of 
(6.1) is closed under pointwise convergence (see, e.g.. Theorem 5.3.2 in [ 8 ]), 
we conclude that u also solves (6.1). This completes the proof of the first 
part of Theorem 1.3. 

In order to prove the second half of the theorem, we keep assuming that u 
is given by (1.5) and we determine the trace of u. For every n, set {Kn, Vn) = 
tr(un). Note that 

Un{x) > uk{x) = -loglPs^iTZnK = 0 ) 

and that uk has trace (A, 0). Indeed, uk is the maximal nonnegative so¬ 
lution of (6.1) that vanishes on dD \ K] see [ 8 ], Theorem 10.1.3. From the 
definition of the trace, it follows that D K. On the other hand, set 


exp(^-y gn{y)X^{dy) 


Ugr.{x) = -logEs^ 

and recall that Ug^ solves (6.1) with boundary condition u\q£) = gn- From 
the bound 




^{nnK= 0 }ew(- J9n{y)X^{dy)] JSn{y)X^{dy) 


<Fs,{nnK ^0) 
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and the previous observations on uk, we see that {un — Ug^){x) converges to 
0 as X ^ y, for every y £ dD \ K. Thus Un has boundary value gn on dD \ K, 
and we conclude that Kn = K and i>n{dx) = gn{x)a{dx). Furthermore, we 
know from Theorem 5.6 of [17] that the convergence of Un to u implies the 
convergence of tr(Mn) towards tr(rt), in the sense of Definition 5.5 of [17], 
and we obtain that tr(ri) = {K,v). 

Finally, if v is any nonnegative solution of (6.1) and {K, p) is its trace, the 
solution u dehned by ( 6 . 1 ) has the same trace as v, and by the uniqueness 
theorem of [17], we must have v = u. 

Remark 6.1. The main contribution of [15] is a direct probabilistic 
proof of the special case /3 = 1 of Theorem 1.3. Note that the probabilistic 
representation of solutions in [15] looks a bit different because it is formu¬ 
lated in terms of excursion measures, which we did not introduce in the 
present work. Very probably (at least in the case d = 2 where the density 
X has a continuous modification), one could give a probabilistic proof of 
Theorem 1.3 along the lines of [15], without any reference to the results of 
[17]. On the other hand, this probabilistic approach remains restricted to 
the values /? < 1, whereas the analytic results hold for any /3 > 0. For this 
reason, we chose to use the full strength of the results of [17] to give a short 
proof of the probabilistic representation (1.5). Also note that closely related 
results appear in the recent work of Dynkin and Kuznetsov; see, for example, 
Theorem 1.4 in [9]. 


APPENDIX 


Proof of Lemma 4.1. First we will prove part (a) of the lemma. From 
the explicit formula (4.1) for the Poisson kernel, we have for every x G L>o 
and yi,y 2 G 

(A.l) lFo(x,yi) - Fo(x,y 2 )l = 7^(1 - \x\^) 

2 vr \yi - x\-^\y 2 - x\^ 

where u ■ v stands for the usual scalar product in M^. Clearly, 1 — jxj^ < 
2 / 9 o(t), and hence, 

\Po{x,yi) - Po{x,y 2 )\P <cpo{x)P 

\yi - x\^P\y 2 - x\^P 


Ei = {xGDo: \yi - x] V 1^2 -x\> 3\yi - ^ 2 ]}, 
E 2 = {xGDo: \yi - x] V 1^2 - x] < 3lyi - y 2 |}- 


Set 
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If X G , we have plainly 

(A.2) |yi-a:| A 1^2 - a;| > 2|yi-y2 |, 

(A.3) \yi - x| A \y 2 - x| > |(|yi - x| V \y 2 - x|). 

Also note that 

\x ■ {yi -y 2 )\ = \ {x - ^{yi + y 2 )) ■ {yi -y 2 )\ 

(A.4) = 1 |(x - yi) • (yi - y 2 ) + (x - y 2 ) ■ (yi - y 2 )| 

< (|x- yi| V |x -y 2 |)|yi -y 2 |. 

By combining (A.2)-(A.4), we obtain 


/El 


Po{xT\Poix,yi) - Po{x,y2)\^dx 


<cf po{xr+p 

JEi 


\x- (yi -y2)P 


■ dx 


< c 

< c 

< c 

< < 


|yi -x| 2 p|y 2 -xpP 

po{x)‘"^P\yi -xr^^|y2 -xr^dxYyi -y2|^ 

\x-ypA\x-y2\>2\yi-y2\ / 

-xl'^^’dx^lyi -y 2 |^ 


^-yi\>‘i\yi-y 2 \ 
2 


([ r^+^+Pr-^Pdr)\y,-y2\P 

\J2\yi—y2\/\2 J 


2\yi-y2\/\2 

c\yi - ^2^, 


if 0 < y < 


2 + a 


“I'- 


c|yi-y 2 p+“-^ 


if y > 


2 + a 

’ 

2 + a 


Then consider the integral on ^ 2 - If x G -E 2 ) we have by (A.4), 

\x- {yi -2/2)1 <3|yi -2/2p- 

Also note that |yi — x| V |y 2 — x| > ^|yi — y 2 |. Then it follows that 

/ po(a;)“|-Po(x,yi)-Po(a;,2/2)rda; 

JE2 


<c [ po{x)'^^P{\yi- x\ ^^ + \y2-x\ ‘^P)dx 

J E 2 

<c f poix)‘^^P{\yi-x\A\y2-x\)~^Pdx 

■/\x-yi\V\x-y2\<3\yi-y2\ 








EXIT MEASURE OF SUPER-BROWNIAN MOTION 


29 


<c/ pQ{xY^'^\yi — x\ dx 

J\x-yi\<3\yi-y2\ 

f3\yi-y2\ 

<c r^+'^-Pdr 

Jo 

<c|2/i-y2p+“-^ 

provided that p <2 a. Hence, the result of part (a) of the lemma follows 
by combining bounds on Ei and E 2 . 

The proof of part (b) is easy. Define 6 = dist(H, 9Do), and recall that 
b> 0. Then from (A.l) we obtain, for every yi,y 2 G QDq, 

sup|Po(T,yi) - Po{x,y 2 )\ <csup ,— <cb~‘^\yi -y 2 \, 

x&B xes |yi-a^rlyi-a^r 

and the result follows. □ 
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